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The full sky cosmic microwave background polarization field can be decomposed into 'electric' and 
'magnetic' components. Working in harmonic space we construct magnetic variables that can be 
measured from observations over only a portion of the sky. Our construction is exact for azimuthally 
symmetric patches, but should continue to perform well for arbitrary patches. For isotropic, uncor- 
rected noise the variables have a very simple diagonal noise correlation, and further analysis using 
them should be no harder than analysing the temperature field. We estimate the tensor mode ampli- 
tude that could be detected by the Planck satellite and discuss the sensitivity of future experiments 
in the presence of a contaminating weak lensing signal. 



Observations of fluctuations in the cosmic microwave 
background (CMB) provide a wealth of information 
about primordial inhomogeneities in the universe. One 
of the most interesting questions is whether there was a 
tensor (gravitational wave) component, as predicted by 
simple inflationary models. CMB polarization measure- 
ments offer a unique probe of this signal ||] . 

Polarization of the cosmic microwave sky is produced 
by electron scattering, as photons decouple from the pri- 
mordial plasma. Gravitational waves produce 'magnetic' 
(i.e. curl) and 'electric' (i.e. gradient) polarization com- 
ponents at a comparable level by anisotropic redshift- 
ing of the energy of photons. Magnetic polarization is 
not produced by density perturbations, so detection of 
a magnetic component would provide strong direct evi- 
dence for the presence of a primordial gravitational wave 
(tensor) component. 

The problem is that the primordial polarization is only 
observable over the region of the sky that is not contami- 
nated by emission from our galaxy and other foreground 
sources of polarization. But the electric/magnetic de- 
composition is inherently non-local, and non-unique in 
the presence of boundaries. To do an exact, lossless sep- 
aration on the incomplete sky one would need to know 
the boundary conditions exactly — both the polarization 
and its derivative around the edge of the observed region. 
Clearly this is not feasible with noisy observed data. 

However the hypothesis of a magnetic signal can still 
be tested. One possibility which avoids differentiating 
the observed data is to construct line integrals of the 
polarization as discussed in Refs. The problem 

with these line integrals is that there are an infinite num- 
ber of them, and they are not statistically independent. 
One would therefore prefer a set of magnetic variables 
to which the electric component does not contribute, but 
which arc finite in number and statistically independent 
(at least for idealized noise properties). In this letter we 
shall construct such variables explicitly for circular sky 
patches. A more complete exposition of our method us- 
ing the spin- weight formalism is given in Ref. p|. 



The observable polarization field is described in terms 
of the two Stokes' parameters Q and U with respect to 
a particular choice of axes about each direction on the 
sky. We use spherical polar coordinates, and the Stokes' 
parameters define a symmetric and trace-free (STF) rank 
two linear polarization tensor on the sphere 

V ab = ^[Q{0®d~4>®4>)-U{d®<i> + <j>®G)]. (1) 

A two dimensional STF tensor can be written as a sum 
of 'gradient' and 'curl' parts 



V ab = V (a V 6) PB - e c ( aV 6) V c P B . 



(2) 



where V is the covariant derivative on the sphere, an- 
gle brackets denote the STF part on the enclosed in- 
dices, and round brackets denote symmetrization. The 
underlying scalar fields Pe and Pb describe electric and 
magnetic polarization respectively and are clearly non- 
local functions of the Stokes' parameters. One can de- 
fine scalar quantities which are local in the polarization 
by taking two covariant derivatives to form V a S7 b 'P ab = 
(V 2 + 2)V 2 P E and e b c V c V a V ab = (V 2 + 2)V 2 P B which 
depend only on the electric and magnetic polarization 
respectively. Here we focus on the magnetic component, 
integrating the latter equation in order to construct an 
observable that depends only on the magnetic component 
of V a b, but does not involve its derivatives. We start by 
defining the surface integral 



B' w = -2 dSW*e b c V c V a P ab , 
Js 



(3) 



where W is a window function defined over some patch 
S of the observed portion of the sky. The factor of minus 
two is included to make our definition equivalent to that 
in Ref. || . Integrating by parts we have 

B' w = V2 f dSW^Vab 
Js 

- 2 I dl a {e b a W*V c V cb - e b c V c W*V ab ) , (4) 

J dS 
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where Wsab = V^^ c ( a ^b) V C W is an STF tensor window 
function. 

We choose window functions so that the line integrals 
that appear in the construction of B' w contain no con- 
tribution from the electric polarization. In general this 
requires that W and \7 a W vanish on the boundary dS. 
The surface integral 

B w =V2 ( dSWf*V ab (5) 
Js 

then depends only on the magnetic polarization. 

We expand the window functions in spherical harmon- 
ics as 

Z>2 \m\<l V V ^ I' 

(The square root factor is included for later convenience.) 
We need not include I — and 1 spherical harmonics 
since they do not contribute to the tensor window func- 
tions. In practice, we are only interested in probing scales 
to some particular / max (the magnetic signal from ten- 
sor modes has maximal power for I ps 100 and decreases 
rapidly with Z), so the sum in Eq. (^|) can be truncated 
at some finite Z max - 

The polarization tensor V a b can be expanded over the 
whole sky in terms of STF tensor harmonics 

Vab = —7= ^2 {Elm ^(hn)afc + B lm Y(l m )ab) I CO 
Im 

where Yff m \ ab and Yn m ^ ab are the gradient and curl tensor 
harmonics of opposite parities defined in Ref. ffl. From 
the orthogonality of the spherical harmonics over the full 
sphere it follows that 

B lm = V2 f dSYg&Vab- (8) 

Jin 

In a rotationally-invariant ensemble, the expectation val- 
ues of the harmonic coefficients define the power spec- 
trum: 

(Bl'm'Blm) — Sl'lSm'mCj 33 . (9) 

The form of the harmonic expansion (^) of the window 
function ensures that the tensor window is 

W B ab = Y, W ^ Y ^)ab, (10) 
im 

where the sum is over I > 2 and \m\ < I. Evaluating the 
surface integrals in Eq. (|^) we find 

B w = Y,W? m Bi m , (11) 

Im 



where the pseudo-multipoles Bi rn are obtained by re- 
stricting the integral in Eq. (||) to the region S: 

Blm = ^ (W+{lm){lm)'Bl> m > - *T^-(im)(im)'^'m') ,(12) 
I'm' 

and the coupling matrices are given by 

W+(Zm)(Zm)' = / &SY^ )ab Yy£y (13) 

** s 

W- (l m )( imy = * I *SYfa )ah Y$*,. (14) 
J s 

The matrix W_n m \n m y controls the contamination with 
electric polarization and can always be written as a line 
integral around the boundary of S. Our aim is to con- 
struct window functions Wi m that remove this contam- 
ination. For azimuthally-symmetric patches the cou- 
pling matrices are block diagonal (M / ±(/ m )(z m )' oc S mm i), 
and W-n m \n m y has only two non-zero eigenvalues for 
\m\ > 2. For m — and \m\ = 1 there are zero and one 
non-zero eigenvalues respectively. 

We now give a practical method for constructing a 
non-redundant set of window functions {M 7 /}, where / 
labels the particular window, that achieve exact separa- 
tion for azimuthal patches. The corresponding (cleanly 
separated) magnetic observables will be denoted Bw t - 
Vectors are now denoted by bold Roman font, e.g. Bw 
has components By/ I , and B has components B\ m , and 
matrices are denoted by bold italic font, e.g. W± have 
components W±n w \u m y. We present the method in a 
form that is applicable (though no longer exact) to arbi- 
trary shaped regions S. 

In matrix form we start with the observed data vector 

B = W+B - iWJE. (15) 

The essence of our method is to choose a window function 
W so as to project out the range of W_ corresponding 
to its non-zero eigenvalues. 

We first diagonalize W+ so it can be written W+ = 
U+D+U\. This allows us to identify the linear combina- 
tions C/+B that arc poorly determined by B — those cor- 
responding to the small diagonal elements of D+ . These 
correspond to polarization patterns that are badly sup- 
ported on the observed patch of the sky and need to be 
removed from the analysis to construct a non-redundant 
set of window functions. The distribution of eigenval- 
ues of W+ is bimodal and the exact definition of 'small' 
is not critical ||. We define an operator U + which 
projects onto the eigenvectors of W+ whose eigenvalues 
are close to one. This amounts to removing the appropri- 
ate columns of C/+ to give the column orthogonal matrix 
C7+. The corresponding smaller square diagonal matrix 
is D+, and we have 

J7|B ps £> + C/|B - iU\WJE. (16) 
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We now multiply by D+ 1 / 2 , denned by [D^^ 2 ]^ = 

Z3^ 1/2 C/|B w f)+ /2 t7|B - iD~ 1/2 UlW-~E. (17) 
This step ensures that isotropic noise in the map gives 

- 1/2 ~t ~ 

isotropic noise on the variables D + J7|B, as shown in 
Ref. §. 

We now project out the unwanted E contamination 
by projecting out of the range of D + ' U + W-. The 
singular-value decomposition for a non-square matrix 
takes the form D^^U+W- = UDV^ where D is diag- 
onal matrix, U is unitary and V is a column orthogonal 
rectangular matrix. There are at most two non-zero sin- 
gular values (diagonal elements of D) per to, and the 
corresponding left singular vectors (columns of U) form 
an orthonormal basis for the range of D+'*U\W-. Wc 
can project out of this range by defining U as the matrix 
obtained by removing the columns of U where the corre- 
sponding singular value is non-zero. Thus, choosing the 
window functions 

W* = &D+ 1/3 Ul, (18) 

we obtain the vector of separated magnetic observables 

B w = W*B = U^D+ 1/2 ulB a U^D^ul B. (19) 

For azimuthal patches the separation is exact; the ap- 
proximation sign arises only from our use of W+ ~ 
in simplifying the matrix that premultiplics 

B. 

In the case of idealized, un-correlated isotropic noise on 
the Stokes' parameters the covariance of the noise ABw 
on Bw is simply Q 

AT = (AB^AB^) = a%I, (20) 

where is the reciprocal weight per solid angle on the 
Stokes' parameters. 

We now consider how to test the null hypothesis that 
the magnetic signal is due entirely to noise. If the noise 
is Gaussian the Bw will be Gaussian and the simplest 
thing to do is a \ 2 test by computing x 2 = BjyiV _1 Bw r 
(for isotropic noise this is just \ 2 = B^B^ ■/ 'ij 2 n ) . Whilst 
the CMB magnetic polarization signal is expected to be 
Gaussian, any spurious or unexpected signal may not be. 
One may therefore also wish to do a more sophisticated 
set of statistical tests at this point. 

If the signal is Gaussian and we know the expected 
shape of the power spectrum we can use the statistic || 

, = B^N^SN^Bw - tv(N~ 1 S) 

^ABlyN^SN^SN^Bw - 2tr(N- 1 SN^S) 

(21) 



m = 




FIG. 1: The real space window functions for an azimuthally 
symmetric sky patch with 9 < 10°, evaluated in the frame 
where the signal is diagonal so the leftmost window produces 
the largest signal for that m. 



giving the number of standard deviations of the maxi- 
mum likelihood from zero for small signal to noise. Here 
the theory covariance is 

s = {b w b w ) = u^b x _l 2 ijXc BB ij+b l _l 2 ij , (22) 

where the diagonal magnetic power spectrum matrix is 
given by [C SB ] (/m ) (/ / m /) = Smm'^u'Cj 38 . The v' distri- 
bution is computed by Monte-Carlo simulation in order 
to assign robust probabilities ||. 

For a particular theoretical model we can rotate to the 
frame where the signal matrix is diagonal, giving rotated 
that are fully statistically independent. In Fig. |l| 
we plot the window functions for the B w which give 
the largest contributions to the signal for a typical flat 
scale-invariant ACDM model. The window functions are 
plotted as line segments of length y/ Q w + U w at an- 
gle tan -1 (Uw I Qw) /2 to the 9 direction where the real 
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FIG. 2: The smallest gravitational wave amplitude At (de- 
fined as in Ref. Q) that could be detected at 99 per cent 
confidence with probability 0.95 (solid lines) and 0.5 (dotted 
lines) by a one year survey that maps a circular patch of sky 
of a given radius assuming uniform noise. The left panel is for 
detector sensitivity s — 10/iK-^/sec (about three times better 
than Planck) and no lensing signal, the right panel is for sensi- 
tivities (bottom to top) s 2 = {0, 25, 50, 100, 200, 400} (iK 2 sec 
and a magnetic lensing signal with white spectrum C^ s = 
4.4 x 10" 6 /iK 2 . 



quantities Qw an d Uw are defined so that 

UB^ = [ dS {Q W Q + U W U) . (23) 

Of the current funded experiments, only Planck is 
likely to detect magnetic polarization if the levels are 
as predicted by standard cosmological models. As a 
toy model we consider the 143 and 217 GHz polar- 
ized channels of the Planck High Frequency Instrument. 
We approximate the noise as isotropic and ignore the 
finite beam width. Combining maps from these two 
channels with inverse variance weighting, we find <jn ~ 
6 x 10~ 3 /iK/K, where Q and U are expressed as dimen- 
sionless thermodynamic equivalent temperatures in units 
of the CMB temperature. We apply an azimuthally- 
symmetric galactic cut of 20 degrees either side of the 
equator and use i max = 250. We find that the null-buster 
would give a detection at 99 per cent significance for a 
tensor initial power spectrum amplitude (defined as in 
Ref. 0) A T > 4 x 10~ 10 with probability greater than 
90 per cent. This corresponds to about 1/10 of the large 
scale Ci detected by COBE, and would be generated by 
inflationary models with slow-roll parameter e > 0.01 and 



energy scale V 1 / 4 > 2 x 10 16 GeV at horizon crossing (for 
example the simplest (j> 2 inflation model). 

Considering experiments more sensitive than Planck, 
one must also take into account the magnetic signal pro- 
duced by weak lensing of the much larger electric polar- 
ization ^ We model this lensing signal as an ad- 
ditional constant Gaussian white noise. In Fig. ^| we 
show the minimum value of At that could be detected 
at 99 per cent significance as a function of the radius of 
the observed patch, assuming a flat ACDM model which 
rcionizes at z = 6.5 (as evidenced by recent Keck obser- 
vations). More sensitive surveys require larger patches 
in order to distinguish the tensor signal from the lens- 
ing signal, and there is an optimal survey size when the 
instrument noise is approximately equal to the lensing 
noise. Even with no instrument noise the lensing signal 
places a lower limit on what can be detected, and in this 
limit one should survey as much of the sky as possible. 
Since we have assumed that source subtraction and com- 
ponent separation can be done exactly our results are 
probably optimistic. 

In summary, we have shown how to construct a set of 
statistically independent magnetic polarization observ- 
ables from observations covering only a portion of the 
sky. For large patches these contain most of the avail- 
able information (see Ref. B for a quantitative discus- 
sion). Some information loss is inevitable as one cannot 
do lossless clean separation without taking derivatives of 
the polarization on the boundary. The variables have 
simple noise properties and can be constructed easily us- 
ing our harmonic-based approach that automatically re- 
moves redundancy due to the finite sky coverage. 
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